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Abstract 

An  iterative  method  is  employed  to  obtain  numerical 
solutions  of  a  two  dimensional  mixed  boundary  value  problem 
for  rectangular  elastic  plates.   The  theory  of  plane  stress 
is  used.   This  problem  is  encountered  when  the  edge  of  a 
uniformly  heated  plate  is  bonded  to  a  rigid  body  and  then 
allowed  to  cool.   The  numerical  results  are  compared  with 
two  formal  asymptotic  expansions  of  the  solution.   They 
correspond  to  the  limiting  cases  of  a  thin  strip  and  a  long 
strip.   The  development  of  bovindary  layers  near  the  edges  of 
the  plate  is  observed.   The  niimerically  determined  stresses 
achieve  their  largest  values  at  the  comer  of  the  bonded 
edge. 
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1.   Introduction. 

A  thin  rectangular  elastic  plate  is  heated  uniformly 
and  then  an  edge  is  bonded  to  a  rigid  half -space.   In  this 
paper  we  present  a  numerical  method  to  determine  the  shrinkage 
stresses  which  develop  when  the  plate  cools. 

We  assume  that  the  deformations  of  the  plate  are 
adequately  described  by  the  theory  of  plane  stress.   Thus  we 
must  solve  a  mixed  boundary  value  problem,  which  we  call 
Problem  S,  for  the  two  dimensional  displacement  equations  of 
equilibrium.   It  is  solved  by  a  previously  developed  [1] 
accelerated  iteration  method.    Each  iterate  is  the  solution 
of  a  mixed  boundary  value  problem  for  Poisson's  equation. 
The  iterates  are  numerically  determined  by  approximating 
this  boundary  value  problem  by  difference  equations.   The 
resulting  system  of  linear  algebraic  equations  is  of  quasi- 
tridiagonal  form.   It  is  solved  by  a  direct  factorization 
method  [4].   The  procedure  is  suitable  for  use  with  a  digital 
computer.   It  is  applicable  to  other  shrinkage  problems  with 
more  complicated  geometries. 

The  method  is  applied  to  rectangular  plates  with 
aspect  ratios,  1/8,  1/U,    1/2,    1.   Without  loss  of  generality 
we  assume  that  the  plate  lies  in  the  x,y  plane  and  that  the 
bonded  edge  coincides  with  the  line  x  =  0,  ly|  ^  a«   Results 
of  the  computations  are  discussed  in  Section  4  and  they  are 


Also  see  references  in  [1]  for  previous  applications  of 
the  method.   Other  applications  are  given  in  [2,5] • 
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shown  in  Figures  I-5.   They  are  compared  with  formal 
asymptotic  expansions  of  the  solution  of  Problem  S  as 
a  — >  00  and  as  a  — >  0.   As  a  —>  00  a  boundary  layer  occurs 
at  y  =  +  a  and  as  a  — >  0  a  boundary  layer  occurs  at  x  =  0. 
The  formal  expansions  are  derived  in  the  Appendix. 

The  stresses  on  the  boiondary  are  discontinuous  at 
the  comers  x  =  0,  y  =  +  a .   The  numerically  determined 
stresses  achieve  their  largest  values  at  and  near  this  point. 
Elsewhere  in  the  plate  the  nonnal  stress  in  the  y  direction 
is  usually  the  largest.  An  approximate  solution  of  Problem  S 
has  been  previously  obtained  [5]  by  truncating  a  formal 
infinite  series  representation  of  the  solution.   The  results 
are  in  qualitative  agreement  with  those  presented  in  this 
paper  except  near  x  =  0^  y  =  +  a.   At  this  point  the  series 
predicts  unbounded  stresses.   The  precise  behavior  of  the 
solution  of  Problem  S  at  x  =  0,  y  =  +  a  is  hot  known. 

2.   Formulation. 

The  mid-plane  of  the  plate  occupies  the  rectangle: 
0<x<l,  |yl  <ain  the  x,y  plane.   The  edge  x  =  0^ 
]y|  <_  a   is  bonded  to  the  rigid  half -space  x  j<  0.   The 
remainder  of  the  boundary  is  stress  free.   Without  loss 
of  generality J  we  assume  that  the  coefficient  of  thermal 
expansion  of  the  plate  is  1.    Since  the  plate  is 
uniformly  heated  the  solution  is  symmetric  with 


respect  to  the  line  y  =  0.   Therefore  we  need  only  consider 
the  rectangle  R:0<x<l,  0<y<a.   We  denote  the 
boundary  of  R  by  B. 

The  two-dimensional  displacement  vector  ufx^y)  has 
the  components  u(xjy)  and  v{x,y).      The  vector  displacement 
equation  of  equilibrium  for  plane  stress  is  [6]; 

Au  =  -  cV(V-  u)  ,     c  =    (l+v)/(l-v)  .  (1) 

and  Hooke's    law   is 

\A.  Jt.  SL. 

o^  =   E    (u^+vvy),     Gy  =   E    (Vy+vu^),     a^y  =    (E    /2)  (l-v )  (u^+v^), 

E      =      E/(l-V^)    . 

Here  E  and  v  are  Young's  modulus  and  Poisson's  ratio  and 
<^  {^>y)>    ^   (^jY)  a-nd  a  (x,y)  are  the  stresses.   Subscripts 
on  u  and  v  denote  partial  derivatives  and  subscripts  on  a 
denote  components  of  the  stress. 

Since  the  plate  is  bonded  to  a  rigid  body,  we  have 
the  boundary  conditions 

u=0,  V  =  y   ,  onx=0,    0<_y<^a.      {j>a) 

The  stress  free  conditions  yield 

u     +VV     =u     +v     =0j  on     x  =   l,      0<y<a,       (5b) 

xyyx  '  >         —J  —      '       ^''    ' 

V     +vu     =u     +v     =0.  on     0<x<l.      y  =  a.       (Jc) 

y'xyx  '  ^     —       '      •}  ^       I 


-3- 


The   symmetry  conditions   are 

u=v=0,  onOj<x<^l,        y=0.  (4) 

The  shrinkage  boundary  value  problem. 

Problem  S,  is  to  determine  two  functions  u(x,y)  and  v(x,y) 
with  continuous  second  partial  derivatives  in  R  that  satisfy 
the  differential  equations  (l)  in  R  and  the  boundary 
conditions  (5). 


5.   The  Numerical  Method. 

We  seek  an  approximate  solution  of  Problem  S  by  an 
accelerated  iteration  method  [1].   Thus  starting  from  an 
initial  estimate  u  (x,y)  of  the  solution  we  define  a 
sequence  of  iterates  u  by  the  recursions 

Au'^  =  -  cV(V  .u^"-"-)  ,  for  x,y  in  R  , 

u"=0,v^=yj  onx=0,   0-<y<a, 


-n 

.  n-1        -n 

n-1 

u     = 

-   Vv          ,    v     = 
y       '     X 

-  u          , 

X 

y 

onx  =  l,   Oj<yj<a, 
Vy  =  -  vuJJ"^  >    "^y  =   -   ^^5'^  '      onO<x_<l,  y=a, 
Uy  =  v"^  =  0  ,  on  0  ^  X  <  1,   y  =  0, 

u"^  =  0u"  +  (l-e)u^'^"-'-  ,  for  x.y  in  R+B. 

Here  9   is  the  acceleration  parameter.   When  9  =  1  we  have 
simple  iterations.   We  find  that  for  other  values  of  9   the 
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(5) 


iterations  converge  more  rapidly.   The  vector  u^  is  a 
provisional  iterate. 

Each  iterate  in  (5)  is  determined  by  solving  the 
following  mixed  boundary  value  problems  for  Poisson's  equation. 

Aw  =  f (x^y)  ,  for  x^y  in  R, 

w  =  g(y)  ,  on  X  =  0,  0  ^  y  <  a. 


(6) 


w^  =  h(y)  ,  on  X  =  1,  0  <  y  <  a, 


X 


Wy  =  p(x)  ,  onO<x<l,  y=a, 

and  on  0  j<  X  j<  I5  y  =  0, 

either   vi     =   0  ,        or   w  =  0. 

Here  f,    g,   h  and  p  are  evaluated  from  the  data  and  from 
previous  iterates. 

We  solve  (6)  by  a  finite  difference  method.   Thus 
R  is  covered  by  a  uniform  rectilinear  mesh  with  mesh  width 
5  in  the  x  and  y  directions  and  the  mesh  lines  are  parallel 
to  the  edges.   The  mesh  is  selected  so  that  the  edges  of 
the  rectangle  coincide  with  mesh  lines.   The  mesh  points 
x.^y.  are  defined  by 

x^  =  i&  ,   i  =  0,1,...,!,    y^.  =  j6,   j  =  0,1,..., J,      (7) 

where  I  =  l/^,    J  =  a/5.   The  set  of  mesh  points  in  R  is 

called  R=  and  the  set  of  mesh  points  on  B  is  called  B=. 

At  each  point  of  R5 +  B=  the  solution  of  (6),  w(x.,y.)j 

is  approximated  by  the  mesh  function  w.  ..   It  is  the 

-'■J 
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solution  of  a  system  of  linear  algebraic  equations  obtained 
by  approximating  (6)  by  appropriate  difference  equations. 
Thus,  the  Laplacian  in  (6)  is  approximated  in  R  by  the 
9-point  Laplace  difference  operator  [7 ]  and  the  boundary 
normal  derivatives  in  (6)  are  approximated  by  three  point 
inward  difference  formulae.   The  derivatives  that  are 
required  to  evaluate  f ,  g,   h  and  p  are  approximated  by 
centered  difference  formulae. 

The  coefficient  matrix  of  the  algebraic  equations 
for  w.  .  is  tridiagonal  with  respect  to  matrices.   The  system 
is  solved  by  a  direct  factoring  method  [4].   Some  of  the 
details  of  the  method  are  described  elsewhere  [1]. 

When  the  iterations  have  converged  for  a  fixed  value 
of  a,    the  stresses  are  computed  at  each  mesh  point  from  finite 
difference  approximations  of  (2).   A  numerical  estimate  of  an 
appropriate  value  of  5  was  determined  from  a  sequence  of  test 
calculations  with  successively  smaller  values  of  5. 

The  following  numerical  convergence  criterion  was 
employed 5 


_  ^.n-li  .  .^-8 


lu^.  -  u^:-^!  <  10-«  .  (8) 


In  a  test  calculation  the  right  side  of  (8)  was  decreased. 
Only  insignificant  changes  in  the  answers  occurred.   The 
condition  (8)  is  only  necessary  for  convergence. 

The  number  of  iterations  that  are  required  for 
convergence  varies  with  a,    6  and  the  choice  of  9.   We  select 
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9   so  as  to  minimize  the  number  of  iterations.   For  most  of 
the  calculations  we  used  9   =    .6.      For  a  =  2  with  5  =  1/20, 
797  iterations  were  required  to  satisfy  (8).   This  used 
approximately  I.5  minutes  of  CDC  66OO  computer  time. 

4.   Presentation  and  Discussion  of  the  Results. 

The  iteration  method  was  applied  to  rectangles 

with  a  =  4,  2,    1,    1/2.   The  mesh  width  5  =  1/20  was  employed 

for  a  =  k,    2,    1,    and  5  =  1/40  was  used  for  a  =  1/2.   Results 

of  the  computations   are  presented  in  Figures  1-J>.      Here  we 

use  the  dimensionless  stresses,  S  =  a  /E  ,  S  =  o  /E     and 

S   =  2a     /FE  (1-v)].   In  these  figures  the  thin  strip  and 
xy     xy '  '•   ^    ^  -'  ^  ^ 

long  strip  approximations  that  are  derived  in  the  Appendix 
are  denoted  by  a  =00  and  a  =  0  .   The  thin  strip  approxima- 
tion is  a  formal  asymptotic  expansion  of  the  solution  of 
Problem  S  as  a  ->  00  .    It  can  only  be  valid  if  y  is  bounded 
away  from  a.   Near  y  =  a  there  is  a  boundary  layer  in  which 


5f 

The  computations  were  performed  on  the  CDC  66OO  Computer 
at  the  AEC  Computing  and  Applied  Mathematics  Center  of 
the  Courant  Institute  of  Mathematical  Sciences.   In  all 
calculations  we  used  v  =  .52. 

We  observe  that  the  plane  stress  theory  may  give  an 
inaccurate  approximation  of  the  three  dimensional  stresses 
as  a  -*- 00  .   In  [8]  it  was  shown  that  the  plane  stress  theory 
gives  the  leading  term  away  from  the  edges  in  a  formal 
asymptotic  expansion  of  the  solution  of  the  three  dimensional 
theory  of  elasticity  as  the  thickness  h  of  the  plate 
approaches  zero,  and  the  other  dimensions  of  the  plate  are 
fixed.   If  the  present  problem  is  considered  as  three- 
dimensional,  there  is  a  double  limit:   h  -> 0  and  a  — >oo . 
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the  solution  of  Problem  S  varies  rapidly  from  the  thin  strip 
approximation  to  satisfy  the  boundary  conditions  at  y  =  a. 
The  long  strip  approximation  is  a  formal  asymptotic  expansion 
as  a  -*-  0.   It  can  only  be  valid  if  x  is  bounded  away  from 
zero.   Near  x  =  0  there  is  a  boundary  layer. 

The  variation  with  y  of  the  stresses  at  the  bonded 
edge  X  =  0  is  shown  in  Figures  1.   We  observe  that  they 
approach  the  thin  strip  approximation 

^x  =  ^xy  =  °  ^    ^y  =  ^  "  ^^  =  '^"^^^   '  ^^^ 

as  a  increases  if  y  is  bounded  away  from  a.   The  stresses 
vary  rapidly  from  (9)  as  y  — >  a.   Since  S  and  S   vanish 

y     xy 

on  y  =  a  (see  Jc)  there  is  a  discontinuity  in  the  stresses 

on  the  boundary  at  x  =  0,  y  =  a.   The  rapid  variation  near 

y  =  a  is  further  illustrated  in  Figure  2  where  the  stresses 

on  the  line  x  =  1/2  are  shown  for  a  =  4. 

S  is  the  largest  stress  over  most  of  the  interior 
y         ^ 

of  the  plate.   The  largest  numerically  determined  stresses  are 
the  values  of  S  near  the  corner  x  =  0,  y  =  a.   The  stresses 
determined  by  the  formal  infinite  series  representation  [5] 
are  unbounded  at  this  point.   Of  course,  the  stresses  that 
are  determined  by  the  present  numerical  method  are  bounded. 
For  the  rectangle  with  a  =  2,  numerical  solutions  were  obtained 
with  the  mesh  widths  6  =  1/20,  6  =  1/26  and  5  =  l/lO.   The 
stresses  were  in  close  agreement  except  near  x  =  0,  y  =  a. 


-8- 


The  stresses  at  this  point  are  shown  in  Table  I. 


Table  I 


6 

1/20 

1/26 

5  =  1/40 

Sx 

4.91-^ 

5.462 

6.454 

S 

y 

2.470 

2.645 

2.956 

-^xy 

4.856 

5.207 

5.825 

Thus  the  stresses  at  x  =  0^  y  =  a  increase  substantially 
as  5  decreases. 

In  Figures  5  we  show  the  variation  with  x  of  the 
stresses  on  the  symmetry  line  y  =  0.   As  a  increases  the 
stresses  approach  (9).   As  a  decreases  the  stresses  approach 
the  long  strip  approximation  near  x  =  1.   Even  for  a  =  1/2 
the  formation  of  a  boundary  layer  near  x  =  0  can  be  observed. 
Presumablyj  the  boundary  layer  will  be  more  pronounced  for 
smaller  values  of  a. 
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Appendix.   The  Thin  Strip  and  the  Long  Strip  Approximations. 
We  define  new  variables  t]   and  U(x,t])  by 

Ti  =  y/a   ,  U(x,T])  =  a"  u(x,y)  ,  (lO) 

where  the  vector  U  has  components  U  and  V.   We  assume  that 
the  solution  U  of  Problem  S  has  the  asymptotic  expansion 

00 

U  -511  U'^(x,Ti)a-J  ,  (11) 

as  a  — >  oo .   The  expansion  coefficients  U*^  ^  j  =  0,1,...^  are 
determined  by  Inserting  (lO)  and  (11)  into  (l)  and  (3)  and 
equating  to  zero  the  coefficlaats  of  each  power  of  a.   Thus 
we  find  that  U'^  and  V^  satisfy  the  following  recursive  system 
of  equations: 

^xx  =  -  ^"x;;^  -  (i+^)in^   '  (15) 

U°  =  0,  V°=  ^,  U'^=  V^=  0,  j  >  0,  on  X  =  0,   0  <  Ti  j<  1,  (14) 

^x  =  -"^n"^  '  ^x  =  -^i''^  >                 on  X  =  1,   0  <  Ti  <  1,  (15) 

^"x  =  "^n'^  '  ^x  =  ~^i~'^  '               on  0  <  X  <  1,  ^  =  1,  (16) 

u^  =  V^'  =  0  ,  on  0  <  X  <1,      ^  =^  0,  {Ij) 

where  we  define  U"^  =  0  if  j  <  0  . 

Setting  j  =  0  in  (12)  and  (15)  we  find  that  U°  and  V° 
are  linear  functions  of  x  with  coefficients  that  are  functions 
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of  r\.      These  coefficients  are  evaluated  from  (l4)  and  (15) 
with  J  =  0.   Thus  we  obtain 

U°  ^  0  ,     V°  =  Ti  .  (18) 

The  bo\indary  condition  (17)  with  j  =  0  is  satisfied  by  (l8). 
We  insert  (l8)  in  (12-15)  with  j  =  1  and  obtain 

U-*"  =  -  vx  ,     V^  s  0  .  (19) 

Equation  (17)  with  j  =  1  is  satisfied  by  (19).   By  induction 
using  (12-15)  and  (18)  and  (19)  it  is  easy  to  show  that 

Ijj   =  Y^   =   0   ,  J  >  2  .  (20) 

The  thin  strip  approximation, 

u^  =  -  vx  ,     v^  =  y  ,  (21a) 

which  we  denote  by  superscript  t,  is  obtained  by  inserting 
(18-20)  in  (11)  and  using  (10).   The  stresses  of  the  thin 
strip  approximation  are  obtained  from  {2),    (10)  and  (21a). 
They  are  given  by 

^x  =  ^xy  =  °  ^    ^y  =  ^-^^-  ^^^^^ 

The  thin  strip  approximation  does  not  satisfy  the 
first  boundary  condition  in  (3c)  at  y  =  a.   Hence  (21)  is  not 
a  uniform  asymptotic  expansion  of  the  solution  of  Problem  S 
as  a  — >oo.   It  can  only  be  asymptotic  to  the  solution  if  y  is 
bounded  away  from  a.   Near  y  =  a  the  solution  of  Problem  S  must 
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vary  rapidly  from  (21)  as  a  — >  oo  so  as  to  satisfy  the 
boundary  conditions  at  y  =  a.   Thus  there  is  a  boundary 
layer  at  the  edge  y  =  a  as  a  — >  oo .   An  analysis  of  the 
boundary  layer  shows  that  the  stresses  in  the  boundary 
layer  are  determined  by  solving  plane  stress  problems  for 
semi-infinite  strips.   These  problems  are  as  difficult  to 
solve  as  the  original  Problem  S.   Thus  we  shall  not  consider 
the  boiindary  layer  approximation.   The  rapid  variation  near 
y  =  a  of  the  solution  of  Problem  S  is  seen      in  Figures  1 
and  2. 

The  long  strip  approximation  is  obtained  by 
substituting  the  change  of  variables  (10)  in  Problem  S  and 
assuming  that  the  solution  has  the  asymptotic  expansion 

00      i 

U  -  ^  V  (x,Ti)a''  ,  (22) 

J  =  0 

as  a  — >  0.   The  V"^  must  satisfy  a  recursive  system  of  equations, 

By  analyzing  this  system  in  the  same  way  as  we  did  (12-17)  we 

find  that 

V^'  =  0  ,      j  >  0  .  (23) 

Thus  the  long  strip  approximation  is  u  •>/  0.   It  does  not 
satisfy  the  second  boundary  condition  (5a)  at  x  =  0.   Hence 
it  can  be  asymptotic  to  the  solution  of  Problem  S  only  for 
X  bounded  away  from  zero.   Near  x  =  0  there  is  a  boundary 
layer.   ¥e  shall  not  analyze  the  boundary  layer.   The  devel- 
opment of  the  boundary  layer  can  be  seen  from  the  numerical 


-13- 


results  presented  in  Figure  3.   The  physical  interpretation 
of  the  long  strip  approximation  is  that  for  small  a,    the 
boundary  constraint  at  x  =  0  influences  the  stresses  only 
near  x  =  0.   Thus  most  of  the  plate  is  essentially  free  to 
contract  as  it  is  cooled  and  no  stresses  are  developed. 
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1.25 


S,(o,y) 


4.914 


6.125 


2.435 


6.376 


-.25 


Fig.  la.   The  variation  of  S  on  the  bonded  edge  x  =  0. 

The  thin  strip  approximation  is  denoted  by  a  =  oo 
and  the  long  strip  approximation  by  a  =  0.   The 
numbers  adjacent  to  the  vertical  arrows  are  the 
numerically  determined  stresses  at  the  edge  y  =  a. 
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Fig.    lb.      The   variation  of  S„  on  the 
bonded   edge. 
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-Sxy(0,y) 


Fig.  Ic.   The  variation  of  S   on  the  bonded  edge. 

xy 
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sa,y) 


Pig.  2.      The  variation  of  the  stresses  along  the  width  of 
the  plate  at  x  =  l/2  for  a  =  4. 
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S,(x,o) 


-.251. 


Fig.  3a.   The  variation  of  S  along  the  length  of  the 
plate  at  the  line  of  symmetry  y  =  0. 
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a  =  00 


S/x,o) 


Fig.  3b.   The  variation  of  S  (x,0) 
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Pasadena,  California  91109 

Prof.  Paul  M,  Naghdi 
Div.  of  Applied  Mechanics 
Etcheverry  Hall 
University  of  California 
Berkeley,  California  9U720 

Prof.  J.  Baltrukonis 
Mechanics  Division 
The  Catholic  Univ.  of  America 
Washington,  D.  C.  20017 

Dr.  Y.  Weltsman 
Div.  of  Applied  Mathematics 
Brown  University 
Providence,  R.I.  02912 


Prof.  A.  J.  Durelli 
Mechanics  Division 
The  Catholic  Univ.  of  America 
Washington,  D.  C.  20017 

Prof.  H.  H.  Bleich 
Departmentof  Civil  Engr. 
Columbia  University 
Amsterdam  &  120th  Street 
New  York,  New  York  1002? 

Prof.  R.  D.  Mindlin 
Department  of  Civil  Engr, 
Columbia  University 
S.  W.  Mudd  Building 
New  York,  New  York  10027 

Prof.  B.  A.  Boley 
Department  of  Civil  Engr. 
Columbia  University 
Amsterdam  &  120th  Street 
Nev  York,  New  York  IOO27 

Prof.  F.  L.  DiMaggio 
Department  of  Civil  Engr. 
Columbia  University 
616  Mudd  Building 
New  York,  New  York  10027 

Prof.  A.  M.  Freudenthal 

Dept.  of  Civil  Engr.  &  Engr.  Mech. 

Columbia  University 

New  York,  New  York  10027 

Prof.  William  A.  Nash 
Dept.  of  Engr.  Mechanics 
University  of  Florida        , 
Jainesville,  Florida  32603 

Prof.  B.  Budiansky 

Div.  of  Engr.  &  Applied  Physics 

Pierce  Hall 

Harvard  University 

Cambridge,  Massachusetts  O2138 

Prof.  P.  G.  Hodge 
Department  of  Mechanics 
Illinois  Institute  of  Technology 
Chicago,  Illinois  60616 
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Universities  (cont'd) 

Prof.  H.  T.  Corten 
321  Talbot  Lab. 
University  of  Illinois 
Urbana,  Illinois  61805 

Prof.  W.  J.  Hall 
Department  of  Civil  Engr. 
University  of  Illinois 
Urbana,  Illinois  61803 

Prof.  N.  M.  Newmark 
Dept.  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois  61803 

Dr.  W.  H.  Avery 

Applied  Physics  Laboratory 

Johns  Hopkins  University 

8621  Georgia  Avenue 

Silver  Spring,  Maryland  209IO 

Prof.  J.  B.  Tiedemann 
Dept.  of  Aero.  Engr.  &  Arch. 
University  of  Kansas 
Lawrence,  Kansas  660k^ 

Prof.  S.  Taira 
Department  of  Engineering 
Kyoto  University 
Kyoto,  Japan 

Prof.  E.  Eeissner 
Dept.  of  Mathematics 
Massachusetts  Inst,  of  Tech. 
Cambridge,  Massachusetts  02139 

Library  (Code  038^1) 

U.  S.  Naval  Postgraduate  School 

Monterey,  California  939i)-0 

Dr.  Joseph  Marin 
Prof,  of  Materials  Science 
Dept.  of  Materials  Sc.  &  Chem. 
U.  S.  Naval  Postgraduate  School 
Monterey,  California  939'j-O 

Prof.  E.  L.  Reiss 

Courant  Inst,  of  Math.  Sciences 

New  York  University 

h   Washington  Place 

New  York,  New  York  IOOO3 


Dr.  Francis  Cozzarelli 
Div.  of  Interdisciplinary 

Studies  and  Research 
School  of  Engineering 
State  Univ.  of  N.Y.  at  Buffalo 
Buffalo,  New  York  1U21I+ 

Dr.  George  Herrmann 
The  Technological  Institute 
Northwestern  University 
Evanston,  Illinois  60201 

Director,  Ordnance  Research  Lab. 
The  Pennsylvania  State  University 
P.  0.  Box  30 
State  College,  Pennsylvania  168OI 

Prof.  Eugen  J.  Skudrzyk 

Department  of  Physics 

Ordnance  Research  Lab. 

The  Pennsylvania  State  University 

P.  0.  Box  30 

State  College,  Pennsylvania  168OI 

Dean  Oscar  Baguio 
College  of  Engineering 
University  of  the  Philippines 
Quezon  City,  Philippines 


Prof.  J.  Kempner 

Dept.  of  Aero.  Engr.  &  Applied  Mech. 

Polytechnic  Institute  of  Brooklyn 

333  Jay  Street 

Brooklyn,  New  York  11201 

Prof.  J.  Klosner 

RSytechnic  Institute  of  Brooklyn 

333  Jay  Street 

Brooklyn, 'New  York  11201 

Prof.  F.  R.  Eirich 

Polytechnic  Institute  of  Brooklyn 

333  Jay  Street 

Brooklyn,  New  York  11201 

Prof.  A.  C.  Eringen 

School  of  Aero.,  Astro.  &  Engr.  Sc. 

Purdue  University 

Lafayette,  Indiana  ^4^790? 
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Universities  (cont'd.) 


Industry  and  Research  Institutes 


Dr.  S.  L.  Koh 

School  of  Aero.,  Astro.  &  Engr.  Sc. 

Puipdue  University 

Lafayette,  Indiana  ii7907 

Prof.  D.  Schapery 
Purdue  University 
Lafayette,  Indiana  U7907 

Prof.  E.  H.  Lee     ^ 
Div.  of  Engr.  Mechanics 
Stanford  University 
Stanford,  California  9^^305 


Mr.  K.  W.  Bills,  Jr. 
Dept.  4722,  Bldg.  0525 
Aerojet-General  Corporation 
P.  0.  Box  19^47 
Sacrarnento,  California  95809 

Dr.  James  H.  Wiegand 

Senior  Dept.  U720,  Bldg.  O525 

Ballistics  &  Mech.  Properties  Lab. 

Aerojet-General  Corporation 

P.  0.  Bojc  19*+7 

Sacramento,  California  95809 


Dr.  Nicholas  J.  Hoff 
Dept.  of  Aero.  &  Astro. 
Stanford  University 
Stanford,  California  914-305 


Prof.  J.  N.  Goodier 
Div.  of  Engr.  Mechanics 
Stanford  University 
Stanford,  California  9^305 


Mr.  J.  S.  Wise 

Aerospace  Corporation 

P.  0.  Box  1308 

San  Bernardino,  California  92U02 


Prof.  Markus  Reiner 

Technlon  R&D  Foundation,  -Ltd. 

Haifa,  Israel 

Prof.  Tsuyoshi  Hayashi 
Department  of  Aeronautics 
Faculty  of  Engineering 
University  of  Tokyo 
BUWKYO-KU 
Tokyo,  Japan 

Prof.  R.  J.  H.  Bollard 
Chairman,  Aeronautical  Engr.  _Dept. 
207  Guggenheim  Hall 
University  of  Washington 
Seattle,  Washington  98IO5 

Prof.  Albert  S.  Kobayashi 
Dept.  of  Mechanical  Engr. 
University  of  Washington 
Seattle,  Washington  98IO5 

Officer- in- Charge 

Post  Graduate  School  for  Naval  Off. 
Webb  Institute  of  Naval  Arch. 
Crescent  Beach  Road,  Glen  Cove 
Long  Island,  New  York  II5U2 


Dr.  Vito  Salerno 

Applied  Technology  Assoc,  Inc. 

29  Church  Street 

Ramsey,  New  Jersey  07^6 

Library  Services  Department 
Report  Section,  Bldg.  lU-lU 
Argonne  National  Laboratory 
9700  S.  Cass  Avenue 
Argonne,  Illinois  6014^0 

Dr.  E.  M.  Kerwin 

Bolt,  Beranek,  &  Newman,  Inc. 

50  Moulton  Street 

Cambridge,  Massachusetts  O2138 

Dr.  M.  C.  Junger 
Cambridge  Acoustical  Associates 
129  Mount  Auburn  Street 
Cambridge,  Massachusetts  02138 

Dr.  F.  R.  Schwarzl 
Central  Laboratory  T.N.O. 
134  Julianalaan 
Delft,  Holland 
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Inci\astry  &  Research  Inst,    (cont'd.) 

Mr.   Ronald  D.   Brown 
jiprpTiied  Physics  Laboratory 
Chemical  Proprulaion  Agency 
8621  Georgia  Avanxie 
Silver  Spring,  Maryland  20910 

Research  and  Development 
Electric  Boat  Division 
General  DynartLcs  Corporation 
Groton,  Connecticut  063iiO 

Supervisor  of  Shipbuilding,   USN, 

and  Naval  Insp.   of  Ordnance 
Electric  Boat  Division 
General  DyDamLcs  Corporation 
Groton,  Connecticut  06^h0 

Dr.   L.  H.   Chen 
Basic  Engineering 
■EleCtiic  Boat  Division 
General  Dynamics  Corporation 
Groton,   Connecticut  O63UO 

Mr.  Ross  H.  Petty 

Technical  Librarian 

Allegany  Balliatics  Lab. 

■Hercules  Powder  Company 

P.  0.  Box  210-  • 

Cumberland,  Maryland  2l501 

Dr.  J.  H.  Thacher 
Allegany  Ballistic .Laboratory 
HercxQes  Powdssr  Company 
Cumberland,  Maryland  21501 

Dr.  Joshua  E.  Green spon 

■J.  G.  Engr.  Research  Associates 

3831  Menlo  Drive 

Baltimore,  Maryland  21215 

Mr.  R.  F.  Landel 
Jet  Propulsion  Laboratory 
U8OO  Oak  Grove  Drive 
Pasadena,  Calif onaia  91103 

Mr-  G.  Lewis 

Jet  Propulsion  'Laboratory 
U8OO  Oak -Grove  Drive 
"Pasadena,  California  91103 


Library 

Newport  News  Shipbuilding  & 

Dry  Dock.  JCon^iany 
Newpoi*t  News,  Virginia  23607 

Mr.  E.  A.  Alex3nder( Research  Dept.(99l) 

Rocketdyne  Division 

North  American  Aviation,   Inc. 

6633  Canoga  Avenue 

Canoga  Park,   California  9130li 

Mr.   Cezar  P.   Nuguid 

Philippine  Atomic  EnSrgy  Commission 
Manila ,   Philippines 

Mr.   S.    C.   Britton 
Solid  Rocket  Division 
Rocketdyne 
P.   0.   Box  51*8 
McGregor,   Texas.  76657 

Dr.  A.    J.    Igna.tQwsJfi 
Redstone  Arsenal  Raffsarch  Div. 
Rohm  &  Haas.  Company 
Huntsville,  Alabama  35807      ' 

Dr.  M.  L.  Merritt 

Division  5412 

Sandia  Corporation 

Sandia  Base 

Albuquerque^   New  Mefxico  87115 

Director 

Ship  Research.  Institute 

Ministry  of  Transportation 

700,  SHINKAWA 

Mitaka 

Tokyo,    JAPAN 

Dr.   H.   N.   Abiaaison- 
Soifthwest  Research  Institute 
8500  Culehira  Road 
San  Antonio,  Teacas  78206 
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Tiidusrtry  &  Renearch  Inst,    (cont'd.) 

Df.   R.   C..   DeHart 

S<Wthwest  Research  Institute 

8500  Culebra  Road 

'San  Arttonio,  Teacas  78206 

Dr.   Thor  Smd.th 

Stanford  Research  Institute 

>Ienlo  Park,    California '  91^025 

Dr.  M.  L-  Baron 

Paul  Weidlinger,  Consulting  Engr. 
777  Third  Ave.  -  22nd  Floor 
New  York,  New  York  10017 

Dr.  T.  C.  Fan 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,  California  904o6 

Dr.  G.  R.  Makepeace 

Director,  Research  &  Engineering 

Lockheed  Propulsion  Company 

P.  0.  Box  ill 

Redlands,  California  9237^ 
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